Melting a stretched DNA 
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We study the melting of a double stranded DNA in the presence of stretching forces, via 3D 
Monte-Carlo simulations, exactly solvable models and heuristic arguments. The resulting force- 
temperature phase diagram is dramatically different for the cases where the force is applied to only 
one strand or to both. Different assumptions on the monomer size of single and double stranded 
DNA lead to opposite conclusions as to whether DNA melts or not as it overstretches. 
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Single molecule experiments have by now enriched our 
knowledge of biopolymer physics at the nanoscale[]J,l3,l3|. 
DNA may be grabbed, twisted and pulled with atomic 
force microscopes, laser tweezers etc., and in this way 
a number of its physical, chemical and elastic proper- 
ties may be measured. Arguably the best characterised 
single molecule experiment is the adiabatic stretching of 
single stranded DNA (ssDNA) and of double stranded 
DNA (dsDNA), the double helical B-form. In most cases, 
the resulting force-elongation curves may be described by 
simple statistical mechanics models, e.g. the worm-like 
or the freely-jointed chain. These inextensible polymer 
models work at relatively low force, and for chain length 
larger than the persistence length, ~ 50 nm (150 base- 
pairs (bp)) for dsDNA and -- 2 nm (4 bp) for ssDNA. 
Still, this classic experiments has prompted several other 
questions, some still outstanding. Most relevant for us 
is the fact that, at about ~ 65 pN, dsDNA elongates to 
about 1.7 times the B-DNA contour length: this regime 
is usually referred to as overstretching. Is overstretched 
DNA melted, ssDNA, or is it a different form of dsDNA, 
named S-DNA [3|, lj|? These two competing pictures are 
not easily discriminated, since the contour lengths of S- 
DNA and ssDNA at the overstretching transition are sim- 
ilar (0.56-0.58 nm per bp). A further complication is 
given by the possibility of pulling one or both strands - 
experiments with laser tweezers may often not be able 
to tell what is happening [3|, as one strand may fall off 
the other due to the presence of nicks breaking the DNA 
backbone, unpeeling [5|. Are the two ensembles equiva- 
lent, or any similar? 

Theories of DNA overstretching are often phenomeno- 
logical, determining parameters from fitting to data [5J, 
lO, l7|. On the other hand, a lot of theoretical work fo- 
cused on simpler models to find the temperature-force 
(T — /) phase diagram of biopoly mers perturbed by ex- 
ternal forces - see e.g. [8|, |9|, [lOl [ill, [l^ ■ DNA overstretch- 
ing was only recently tackled within similarly simplified 
models, by assuming that it amounts to melting ds to 
ssDNA [13|, llJ] . On the basis of a continuum model [l3| , 
it was proposed that when one strand is being pulled. 



there should be a "reentrance" at low T, similar to that 
found for DNA unzipping [9|. Another interesting de- 
bated question is whether the ds phase is destabilised or 
not by a stretching force [3|, ll3|. A final prediction is 
that the melting transition, which is first order at / = 0, 
should become second order at f =^ 13lll4|. 

Here we will be concerned with the melting transition 
of a model of dsDNA anchored on one side and subject to 
external stretching forces, fi and £2 on the two strands on 
the other side (Fig. 1). At variance with previous studies, 
we introduce two distinct pulling modes: (A) f 1 = f and 
£2 = 0, i.e. one of the two strands is stretched and the 
other one is left free, DNA unpeeling; (B) f 1 = f2 = f/2, 
i.e. both strands are stretched by independent, but equal, 
forces, stretched DNA denaturation. Each strand will 
be considered at some coarse grained level and modeled 
as a self-avoiding polymer. Strand complementarity is 
taken into account by allowing a pair of monomers in 
different strands to be bound with an energy — e < only 
if their monomer index along the two strands is equal [l5| . 
Monomer sizes (length of 1 bp) and persistence lengths 
are here the only parameters controlling the microscopic 
properties of dsDNA and ssDNA. We choose monomer 
sizes in two fundamentally distinct ways. If we assume 
dsDNA exists also in the stretched S-form, monomer sizes 
are then equal for both ssDNA and dsDNA (cases Aq Bq). 
If not, dsDNA monomer size, from the B-form, is smaller 
than for ssDNA (cases Ai Bi). The values of monomer 
sizes and persistence lengths mostly impact on the large 
and small force behaviour of the system respectively. 

From now on T will be measured in units of fcs/e {ks 
is the Boltzmann constant). Unlike all previous work, 
we perform numerical simulations finding a phase dia- 
gram valid in the whole T range. We can thus discover 
that models A and B behave amazingly differently, both 
qualitatively and quantitatively. 

We begin with heuristic arguments to determine the 
critical force, fc{T), separating the zipped phase from 
the unzipped/melted one. At / = we have the standard 
melting phase transition occurring [l^, [l6| at T = Tm, 
i.e. fc{Tm) = 0. For case Aq, unpeeling with same 




FIG. 1; Schematics of our model. A dsDNA molecule is an- 
chored at one of its ends, whereas two pulling forces stretch 
each strand at the other ends. At a finite T bubbles may 
appear in typical DNA configurations as in the figure. 
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monomer size for ss and dsDNA, at large / the (almost) 
eompletely stretehed strand acts as a ID substrate for 
the adsorption of the free (un-stretched) strand. In this 
case the only allowed conformation are Y-like configura- 
tion with no bubbles. An energy-entropy argument leads 
to a vertical asymptote for fciT) at Ta = e/ss, where s.. 



is the entropy per monomer of a single DNA strand [17 1 . 
This argument is fully confirmed by numerical results 
coming from 3D Monte-Carlo simulation of model (A), 
i.e. two paired DNA strands, one of which is under a 
stretching force of modulus / (see Fig. 2, with typical 
snapshots from the simulations shown as well). Each 
strand is modelled by a self-avoiding chain, made up by 
N beads of size ct (~ 1 nm, the size of ssDNA), con- 
nected by springs. To model homologous base pairing, 
we chose a truncated Lennard-Jones potential, with min- 
imum value — e, attained when the corresponding beads 
are at a mutual separation a. The bonds between succes- 
sive beads are harmonic springs with minimal and max- 
imal elongations equal to 0.7 a and a respectively (the 
maximum elongation is then {N — l)a). The simulations 
were performed by proposing local deformations of the 
chain. We used multiple Markov chains [18[ to improve 
sampling efficiency at small T or large /, and reweighted 
the simulation data prior to analysis according to the 
Ferrenbcrg-Swcdscn algorithm [l9|. The critical points 
were estimated via the peaks of the specific heat. Our 
data suggest a first order transition at all / ^ 0. 

The phase diagram found numerically (Fig. 2) shows 
that there is as expected a vertical asymptote. At small 
forces, the unpeeling pulling mode stabilises the zipped 
dsDNA, so that the critical temperature, Tc{f), increases 
with / (this agrees with a thermodynamic analysis of 
DNA overstretching data [i,[il). 

To study model Bq, we have performed our simula- 
tions in the case in which both DNA strands are under 
the action of the same force (/i = /2 = //2). As bubbles 
(see Fig. 1) are stretched out, the two strands are on 
average closer to each other with respect to the f — 
case and one may argue that the dsDNA should be sta- 
bilised by /. Another expectation is that at very large 
/, both strands become straight in the same direction 
and the zipped phase is again stabilised since base pair- 
ing does not lead to further entropy cost. Both these 
expectations are confirmed by the Monte-Carlo simula- 



FIG. 2: Solid line: f—T phase diagram in the continuum limit 
for a model with two self-avoiding chains kept together by a 
truncated Lennard-Jones potential (with minimum — e), and 
a force stretching one of the strands. There is a first order 
transition separating the zipped from the unzipped phases. 
Dashed line: same for a force stretching both strands. The 
temperature range includes unphysical values to show the 
whole theoretical diagram. 



tions (see Fig. 2, dotted line): the region of stability 
of the zipped phase, in the (T, /) plane, increases, and 
the denaturation temperature appears to be a monoton- 
ically increasing function of /. The phase diagrams for 
both models Aq and Bq are somewhat surprising as at 
low T (and indeed at any T for model Bq) DNA can 
stay zipped even at large /. This would however be com- 
patible with the interpretation of overstretched DNA as 
a ds form different from B-DNA [J]. The key point is 
that monomer sizes of S-DNA and ssDNA are roughly 
the same. The overstretching transition (B-DNA to S- 
DNA) takes place in Fig. 2 at 65pN (the energy scale 
e was set to 2.7 Kcal/mol, by assuming T„i = 70 °C). 
For unpeeling, overstretching can be followed by further 
unzipping to ssDNA, whereas this is not possible when 
pulling both strands. A second elongation transition is 
indeed observed in some cases [5|, |7[. Our model also 
suggests that at higher T, and when pulling just one 
strands, overstretching should change to a melting transi- 
tion, even assuming the existence of S-DNA. The picture 
emerging from our results might thus reconcile both in- 
terpretations of overstretched DNA. 

It is instructive to compare the phase diagrams found 
with the Monte-Carlo simulations to those of an exactly 
solvable model for DNA melting under a stretching force. 
In these simplified models the two strands are mutually 
avoiding directed walks along the [11] diagonal of a 2D 
lattice (Fig. 3a). When two beads arc in contact they 
gain a pairing potential e, and one end of the chain is 
fixed at the origin for both strands, while at the other 
end there are stretching forces, /i and /2 (Fig. 3b). The 
case with /2 = corresponds to our model Aq. 

One can write down a recursion for the partition func- 
tion of a system with only one degree of freedom, the 
open end distance x. If the strands are subjected to two 
forces /i and /2 (e.g. identified by their projection along 




FIG. 3: (a) Geometry for the exact calculations, (b) Phase 
diagram for unpeeling (solid line) and stretched denaturation 
(dashed line), (c) Fraction of zipped bases as a function of T 
for / = (solid line), for DNA unpeeling, and / = 1 (dashed 
line) , for stretched denaturation and / = 2 (dot-dashed line) . 
(d) End-to-end distance of the pulled strand per monomer 
as a function of / for DNA unpeeling and T — 2 (solid line), 
T — 3 (dashed line) and for stretched denaturation and T — 4 
(dot-dashed line, with the arrow denoting the transition). 



the positive [-1 1] vector, Fig. 2a), then these are: 

Zn+i^x) = Zn^x + I)y2/yi + ZNix - l)yi/y2il) 
+ Ziq{x){yiy2 + ll{yiy2)) forx > 

ZN+me'^'^ = ZN{l)y2/yi (2) 

+ pZN{Q){yiy2 + l/{yiy2))- 

In Eq.s lll2l ZAf(a::) indicates the partition function of two 



walks with open end distance x, yi — €■ 
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e^^f'^ , and 



p is a parameter which controls the weight of dsDNA seg- 
ments and which we temporarily set to 1, implying that 
both monomer sizes and persistence lengths are equal in 
both ss and ds phases. The phase diagrams for models 
Af) and Bq can be found with e.g. a generalisation of the 
methods used in Ref. Q and are plotted in Fig. 3b. 

Fig. 3b shows that the phase diagram is qualitatively 
similar to the one already found in 3D Monte Carlo sim- 
ulations in the continuum. The presence of a vertical 
asymptote for unpeeling, &i T = T^ = r^, the tran- 
sition temperature for polymer adsorption on a wall, is 
confirmed. Interestingly, while unpeeling, at / 7^ is first 
order, stretched denaturation is second order, as can be 
seen from the plots of average percentage of zipped bases 
(0, Fig. 3c), and of the average open end distance per 
monomer ((a;)/N, Fig. 3d). (The unpeeling phase dia- 
gram has been found, in another context, in [12l|.) Sim- 
ulation data show consistently a much smoother transi- 
tion when pulling both strands. On the other hand, in 
the case of DNA unpeeling there are some differences be- 
tween Monte Carlo simulations and exact results, most 
notably the shape of the phase diagram close to r,„ , as 
the critical temperature decreases with f. 

Interestingly, an expected difference concerns the order 



of the melting transition at zero force, which is thought 
to be first order in 3D [16[, and is second order in the 
directed model. Another realistic feature that the simu- 
lated model reproduces is the larger effective persistence 
length of the ds chain with respect to the ss one, at vari- 
ance with the lattice model. 

To gain more insight, we present a Flory-like argument 
to rationalise the / — > behaviour of the phase diagram 
for DNA unpeeling. We consider both cases ^0 and At, 
in which the effective length of dsDNA segments [20| 
is larger (Ed - 300 bp ~ 100 nm for B-DNA) with re- 
spect to that of ssDNA ones (€s '^ 8 bp ~ 4 nm). The 
monomer sizes are ad for dsDNA and as — 0.56 nm for 
ssDNA, and may differ. Let us consider T ~ Tm and a 
small force / ^ 1 . If we disregard denaturation bubbles 
(considering only Y-shaped configurations), the energy 
of two DNA strands with m = ON open bp, under the 
action of an unpeeling force f , is i? = —e{N — tti) — f • r, 
where r is the position of the un-anchored end of the 
stretched strand. We may expect the probability den- 
sity, F(r|iV, m) for the stretched strand to be described 
by a Gaussian approximation, in _D-dimension, 



P{r\N,m) = (27rcr2(6))iV 
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The total number of Y-like configurations is given by, 
modulo power law corrections, C{N\m) oc exp{ ^ Z" '^"^ 
jig} where the entropies per monomer of the double 
and single stranded DNA have been estimated to be 
s/id and s/is respectively. The total entropy of Y- 
like configuration with a fixed r is estimated as 5 = 
\og{C{N\ra)P{Y\N,ra)}. Minimizing the free energy 
T ^ E — TS with respect to r and 6, we obtain the 
following form for the critical force, fc{T), 
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Thus iiiddd > ^s«s, as in reality, Tc increases with /. No- 
tice that fdT) ~ (T - r„) 1/2 at T ~ r+ (self-avoidance 
would change the exponent from 1/2 to ;/ ~ 0.588, al- 
though this trend is not discernible in Fig. 2 due, possi- 
bly, to the distance from the thermodynamic limit). The 
behavior of fc given by eq.(l5]) is only valid for T '^ Tm- 
Away from T,„ this transition line has to turn left and 
e.g. approach a vertical asymtote for Og = ad- This 
is indeed the shape of the phase diagram obtained with 
Monte-Carlo simulations for model Aq with a,, ~ ad, as 
£d > is due to the interaction. Why is the dsDNA phase 
instead destabilised by / in the exact model (fig. 3b)? A 
plausible reason is that the argument we have just pro- 
posed neglects denaturation bubbles, and thus predicts 
T„j = Ta for all forces, when £d = is, o,d = o-s- Bub- 
bles are unimportant for the continuum simulations, as 
the melting is first order for self- avoiding chains, but are 



crucial for the exact model, as the transition is second 
order and there are ~ N^''^ bubbles at the transition [9|. 
Bubbles increase T„ and the stability of the zipped phase 
at / = 0, and might indeed cause the small force reen- 
trance to disappear. 

We now turn to models Ai and _Bi, assuming that 
dsDNA can only be in the B-form. It is then possible 
to generalize our exact results by introducing different 
values ^ ~ ^d/^s ~ 25 and a = aa/O's = 0.6 for Kuhn 
lengths and monomer sizes of B-DNA and ssDNA. To in- 
corporate these elements into our models, we may weigh 
dsDNA segments in Eq. [2] as follows: 

p(/3,/,e,a) - [2cosh(/3M)]'/V2cosh(/3/) (6) 




FIG. 4: (T, /) phase diagrams for model A-i (solid line) and 
B\ (dashed line), with ad — 0.34, Id = 50 nm, and a^ — 0.56, 
^3 = 2 nm for ds and ss DNA respectively. 



The resulting phase diagrams are shown in Fig. 4. The 
most striking difference with respect to the phase dia- 
grams in Figs. 2 and 3 is that there is now a bound 
region for the zipped phase, and at large force the DNA 
is always ss and denatured, in agreement with the in- 
terpretation of overstretched DNA as melted DNA. A 
second difference is that close to the denaturation tran- 
sition the shape of the critical line is affected by denatu- 
ration bubbles as well. As a result for directed walks ds- 
DNA is stabilised by a force only when both strands are 
stretched. When the transition is first order (DNA un- 
peeling), then we can derive dfc{T)/dT = —SSu,z/Sxu^z, 
where SSu,z and dxu,z are the changes in entropy and 
elongation (along the direction of f ) between the zipped 
and the unzipped phases [31] . Close to T = for DNA un- 
peeling both these quantities are positive, hence there is 
no reentrance, whereas close to the melting temperature 
SSu,z > but the sign of Sxu^z niay vary. 

To conclude, we presented Monte-Carlo simulations 
and exact results to determine the phase diagram of 
dsDNA melting in the presence of an external force, 
stretching one or both strands ~ DNA unpeeling or 
stretched DNA denaturation. Our results show that 
these two cases, which may often be difficult to distin- 
guish in single molecule experiments, lead to strikingly 
different results, both qualitatively and quantitatively. 
Contrarily to previous claims, our results suggest that 
DNA unpeeling is a first order phase transition. We have 
discussed how to introduce in our models the realistic 
persistence lengths and monomer sizes of ds and ssDNA. 
This can be done in different ways according to what 
is assumed about the nature of the overstretched DNA 
state. Even if we assume dsDNA can exist in the S-DNA 
form we have shown that overstretching transition to 
a melted ssDNA is possible when one strand is being 
pulled. According to our calculations, we may infer the 
phase diagram of a realistic model, taking these effects 
as well as self-avoidance, fully into account. We predict 
that it should (i) display no "reentrance" for low T, and 
(ii) show dsDNA stabilisation by force, close to Tm 21 1. 
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